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Self-diffusion of a polymer chain in a melt is studied by Monte Carlo simulations of the bond 
fluctuation model, where only the excluded volume interaction is taken into account. Polymer 
chains, each of which consists of N segments, are located on an L x L x L simple cubic lattice 
under periodic boundary conditions, where each segment occupies 2x2x2 unit cells. The 
results for N = 32,48,64,96,128,192,256,384 and 512 at the volume fraction <j> ~ 0.5 are 
reported, where L = 128 for TV < 256 and L = 192 for TV > 384. The TV-dependence of the self- 
diffusion constant D is examined. Here, D is estimated from the mean square displacements 
of the center of mass of a single polymer chain at the times larger than the longest relaxation 
time. From the data for TV — 256, 384 and 512, the apparent exponent Xd, which describes the 
apparent power law dependence of D on TV as D oc N~ Xd , is estimated as %d — 2.4. The ratio 
Dr/{Re) seems to be a constant for TV = 192, 256, 384 and 512, where r and (B%) denote the 
longest relaxation time and the mean square end-to-end distance, respectively. 
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Recently, many works have been done in order to study 
the slow dynamics of a single polymer chain in a melt 
through simulations. 1_9 ' These studies are motivated 
by the reptation theory of the dynamics of concentrated 
polymer systems. 10 ' 11 ) In the reptation theory, the dy- 
namics of concentrated polymer systems is described in 
terms of the motion of a polymer chain trapped in a 
tube, which represents the entanglement effects of the 
surrounding polymers. The self-diffusion constant D of 
' the center of mass and the longest relaxation time r of a 
polymer chain of TV segments are predicted to behave as 
D oc TV~ 2 and t oc TV 3 , respectively. On the other hand, 
, the experimental results on the behaviors of D and r 
' have been summarized as D oc TV~ 2 and r oc TV 3 - 4 . 10 ' 11 ' 
Thus, the exponent for r observed in the experiments 
is different from that predicted by the reptation theory. 
, This discrepancy between the experiments and the repta- 
' tion theory is considered to be attributed to the contour 
length fluctuation. 10, 12 ~ 14 ' The contour length fluctua- 
' tion makes the relaxation time r shorter than that pre- 
, dieted by the original reptation theory. As N/N c be- 
comes very large, the effect of the contour length fluctu- 
ation becomes negligible and r approaches the predicted 
power law behavior TV 3 from below, where TV e denotes 
the number of segments between entanglement points. 
Because of this crossover behavior, the apparent expo- 
nent for t increases beyond three and reaches 3.4 for 
intermediately large values of N/N e . This behavior of 
the apparent exponent for r has been confirmed by the 
present authors through Monte Carlo simulations. 1 ' 

In contrast to the relaxation time r, the exponent for 
the self-diffusion constant D observed in the experiments 
has been believed to agree with that predicted by the 
reptation theory. Therefore, the contour length fluctu- 
ation has been considered to have little effect on the 
TV-dependence of the self-diffusion constant D. In the 
recent experiments by Tao, Lodge and Meerwall, 15 ' 16 ' 
however, it is reported that the self-diffusion constant D 



in hydrogenated polybutadiene (hPB) alkane solutions 
and hPB melt behaves as D oc TV~ 2 ' 4±0 . Theoretically, 
Frischknecht and Milner 17 ' reported that the apparent 
power law dependence of D on TV is modified by the con- 
tour length fluctuation as D oc TV" 2 - 25 at TV/TV C ~ 80, 
while it approaches the behavior D oc TV~ 2 predicted 
by the reptation theory as TV/TV becomes larger. Thus, 
the behavior of the self-diffusion constant D is still in 
dispute. Therefore, it is desirable to elucidate the TV- 
dependence of D through simulations. Even in the re- 
cent simulations of a polymer melt, 7 ' 8 ' however, the 
TV-dependence of D has not been clarified. Dunweg, 
Grest and Kremer 7 ' performed molecular dynamics sim- 
ulations for TV = 5, 10, 25, 40, 50, 100, 200 and 350 at the 
reduced density p — 0.85 and the reduced temperature 
k B T = l. 5 ' 6 ' They estimated D only for TV < 200. From 
their values of D for TV = 100 and 200, the apparent 
exponent Xd in the apparent power law D oc N~ Xd can 
be estimated as Xd — 1.93. However, the values TV = 100 
and 200 are not in the range where the apparent expo- 
nent for the relaxation time r exceeds three. 7 ' There- 
fore, their results cannot resolve the dispute as to the 
apparent power law dependence of D on TV. Moreover, 
the sizes of the systems used in their simulations do not 
seem to be sufficiently large. The ratio of the linear di- 
mension of the simulation box to the root mean square 
end-to-end distance of a polymer chain is 1.76 and 1.54 
for TV = 100 and 200, respectively. Kreer, Bashnagel, 
Muller and Binder 8 ' performed Monte Carlo simulations 
of the bond fluctuation model 18 ' for TV = 16, 32, 64, 128 
and 512 at the volume fraction <fi = O.5. 3 ' From their 
values of D for TV = 128 and 512, the apparent expo- 
nent Xd is estimated to be a little smaller than two. Al- 
though the values TV = 128 and 512 are in the range 
where the apparent exponent for the relaxation time ex- 
ceeds three, 1 ' 8 ' they are only two values of TV and it is 
not clear how the apparent exponent Xd changes as TV 
becomes large. Furthermore, the problem of the system 
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size also exists in their simulations. The ratio of the 
linear dimension of the system to the root mean square 
end-to-end distance is 1.75 for TV = 512. 

The purpose of the present letter is to clarify the TV- 
dependence of D in the range of values of TV where the 
apparent exponent for the relaxation time exceeds three. 
As in our previous letter, 1 ' we perform Monte Carlo sim- 
ulations of the same model as used in Ref. 8. More values 
of TV than in Ref. 8 are studied. For each value of TV, 
the linear dimension of the simulation box is at least 2.5 
times larger than the root mean square end-to-end dis- 
tance. The self-diffusion constant D is estimated from 
the mean square displacements of the center of mass of 
a polymer chain at the times larger than the longest re- 
laxation time. 1 ' 

In the present study, the bond fluctuation model 18 ' is 
used for Monte Carlo simulations of a polymer melt. A 
system of M polymer chains, each of which consists of TV 
segments, is considered on a simple cubic lattice of which 
the lattice constant is a. Each side of the lattice has L 
lattice points and the periodic boundary conditions are 
applied. The polymer chains are represented by M series 
of TV lattice points {(ri >m , ■ ■ ■ ,r"jv, m ); m — 1, • • ■ , M}, 
where r^ m denotes the center of the ith segment of the 
mth polymer chain. Each segment occupies a cube of 
2 3 unit cells and no overlap of segments is allowed. The 
volume fraction is given by <f> = 2 3 MN/L 3 . The lengths 
of the bond vectors are restricted to the five prescribed 
lengths {2a, V5a, y6a, 3a, ylOa}. The center of each 
segment stochastically moves to one of its nearest neigh- 
bor lattice points with the transition probability per unit 
time 7/6 if the new position is allowed by the conditions 
for the excluded volume and the bond vectors. The dy- 
namics is simulated by a standard Monte Carlo method, 
which uses discrete time steps and updates MTV seg- 
ments in random sequence. The time scale 7 -1 corre- 
sponds to one update per segment. 19 ' We choose a = 1 
and 7 = 1. 

In the present study, the cases of the vol- 
ume fraction (f> ~ 0.5 are examined for TV = 
32, 48, 64, 96, 128, 192, 256, 384 and 512. We choose L = 
128 for TV < 256 and L = 192 for TV > 384. For these 
values of L, it is confirmed that no polymer chain inter- 
acts with its images generated by the periodic boundary 
conditions in the present simulations. The vectorized 
code for the bond fluctuation model 20 ' is used on vec- 
tor processors for TV > 192. The equilibrium average is 
calculated as the average of the time averages obtained 
from Monte Carlo simulations of four statistically inde- 
pendent systems. In each Monte Carlo simulation, the 
time average is calculated from the sequence of Nj states 
taken at intervals of t\ after the initial period of 2j, which 
is discarded for the equilibration. The values of the pa- 
rameters used in the present simulations are shown in 
Table I. 

Figure 1 shows the elapsed time dependence of the 
mean square displacement ([r c (i) — r c (0)] 2 ) of the cen- 
ter of mass r c of a polymer chain, where t denotes the 
elapsed time. In Fig. 1, the mean square displacement 
and the elapsed time are scaled by the mean square end- 
to-end distance (i? 2 ) and the longest relaxation time r, 
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Fig. 1. The elapsed time t dependence of the mean square dis- 
placement {[r c (t) — r c (0)] 2 ) of the center of mass of a polymer 
chain. The mean square displacement and the elapsed time are 
scaled by the mean square end-to-end distance (Rg) and the 
longest relaxation time r, respectively. 
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Fig. 2. Log-log plot of DN versus N. Solid line represents the 
result of the fit of the data for N = 256, 384 and 512 explained 
in the text. The results of Ref. 8 are also shown for comparison. 
Dashed line represents the theoretical result D oc N~ 2 . 



respectively. The values of (i? 2 ) and r are shown in Ta- 
ble I. For TV < 256, the values estimated in Ref. 1 are 
used, while those estimated from longer simulations than 
in Ref. 1 are used for TV > 384. Here, the longest relax- 
ation time t is estimated by the same method as in Ref. 
1, which is based on the variational estimation of the re- 
laxation modes and rates 21 ' and has been applied to var- 
ious polymer systems. 1, 19 ' 22 ~ 25 ' The self-diffusion con- 
stant D is estimated by fitting the data points at times 
longer than the longest relaxation time r to a straight 
line ([r c (t) — r c (0)] 2 ) = 6Dt + constant. The estimated 
values of D are shown in Table I. 

Figure 2 shows a log-log plot of DN versus TV. The 
apparent exponent x& of the power law dependence of 
the self-diffusion constant D cx N~ Xd , which appears as 
the slope in the log-log plot, increases as TV increases. 
It is smaller than two for TV < 96 and larger than two 
for N > 128. The least squares fit of the data points for 
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Table I. The parameters of the simulations. The mean square end-to-end distance (R^.), the inverse of the longest relaxation time 1/r 
and the self-diffusion constant D of a polymer chain are also shown. The values of {R^) and 1/r for TV < 256 are taken from Ref. 1. 
The symbols k and M denote xlO 3 and xlO 6 , respectively. 
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N = 128, 192 and 256 to a straight line in a log-log plot of 
D versus N gives x d ~ 2.21, while that for N = 256, 384 
and 512 gives Xd — 2.44. The latter result is shown in 
Fig. 2 as a solid line. 

The reptation theory predicts that r oc iV 3 and D oc 
iV -2 in the large N limit, while (R 2 ) oc N is expected to 
hold in a melt. Thus, the ratio Dt/(R 2 } is predicted to 
be independent of N for sufficiently large N. From the 
values in Table I, the apparent exponents of the power 
law dependences (R 2 ) oc N 2v and r oc N x ' are estimated 
as 2v ~ 1.02 and x T ~ 3.52 for N = 256, 384 and 512. 
As mentioned before, the apparent exponent Xd in D oc 
N- Xd is estimated as x d ~ 2.44 for N = 256, 384 and 
512. Thus, the relation — x d + x r — 2v ~ seems to hold. 
This means that Dt/(R 2 ) is independent of N even in 
the range of N where the true asymptotic behaviors of 
r and D are not expected to be seen. This behavior can 
be seen from the fact that the slopes in Fig. 1, which 
correspond to 6£>r/(i? 2 ), are almost the same for N > 
192. 

According to the reptation theory, the ratio Dt/(R 2 ) 
is not only independent of N but also a constant which 
contains no adjustable parameter. Its value is given by 
l/(37r 2 ) ~ 0.034. 10 ) Figure 3 shows a semi- log plot of 
Dt/(R 2 c ) versus N. For N > 192, the values of Dt / (R 2 C ) 
seem to converge to a constant value around 0.024. If this 
value can be regarded as the large N limit of Dt/(R 2 ), 
the value is about thirty percent smaller than the value 
predicted by the reptation theory. The agreement be- 
tween the two values is considered to be reasonable if we 
take into account the following facts. Firstly, the esti- 
mated value of t is considered to have a tendency to be 
smaller than the exact value, because the relaxation time 
r is estimated by the variational method. 1 ' 21 ' Secondly, 
the accuracy of the predicted value of Dt/(R 2 ) itself is 
unknown, since the reptation theory is based on a simple 
model of polymer melts. 

In summary, we have examined the -ZV-dependence of 
the self-diffusion constant D of a polymer chain of N 
segments in a melt through Monte Carlo simulations of 
the bond fluctuation model. The apparent exponent Xd 
in the power law dependence D oc N~ Xd is found to 
increase beyond two as N increases. From the data for 
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Fig. 3. Semi-log plot of Dt/{R%) versus N. Solid line represents 
the average of the values of Dt/(RI) for N = 192, 256, 384 and 
512. The results of Ref. 8 are also shown for comparison. In the 
results of Ref. 8, r v — i 8 ) is used as the longest relaxation time. 



N = 256, 384 and 512, x d is estimated as x d ~ 2.44. 
For large values of N, the ratio Dt/(R 2 ) seems to be 
a constant value around 0.024, which reasonably agrees 
with the theoretical prediction. The reptation theory 
predicts that Dt/ (R 2 ) is independent of the monomer 
concentration of a melt. 10 ) Therefore, it is interesting to 
investigate the volume fraction <fr dependences of D, t 
and (i? 2 ) in addition to their iV-dependences. A study 
in this direction is in progress 
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